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Abstract 



In classical continuum mechanics, quasi-linear systems of conservation 
laws can be symmetrized if they admit an additional convex conservation 
law. In particular, this implies the hyperbolicity of governing equations. 
For capillary fluids, the internal energy depends not only on the density 
but also on its derivatives with respect to space variables. Consequently, 
the governing equations belong to the class of dispersive systems. In that 
case we propose a symmetric form of governing equations which is different 
from the classical Godunov - Friedrichs - Lax representation. This new 
symmetric form implies the stability of constant solutions. 

1 Introduction 

Quasi-linear systems of conservation laws can be symmetrized, if they admit 
an additional convex conservation law (Godunov, 1961, Friedrichs and Lax, 
1971). The symmetric form implies hyperbolicity of governing equations. For 
conservation laws with vanishing right-hand side, the hyperbolicity is equiva- 
lent to stability of constant solutions with respect to perturbations of the form 
e «'( k * x -At) j j2 _ _^ k* _ ^ . . . x _ . . . >x n^ wn ere "*" denotes the 

transposition. Indeed, the following symmetric form of governing equations for 
an unknown vector variable v 



where matrix A = A* is positive definite, B % — (B 1 )* , implies the dispersion 
relation 
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which determines real values of A for any real wave vector k. 

In this note we get an analog of symmetric form (1) for equations of capillary 
fluids that belong to the class of dispersive systems, because the internal energy 
depends not only on the density but also on its derivatives with respect to space 
variables. We will see that the analog of equation (2) is 

det(B + iC - XA) = (3) 

where C = — C* is an antisymmetric matrix depending on the wave vector 
k. Since B + iC is Hermitian matrix and the symmetric matrix A is positive 
definite, all the frequencies A are also real. For a capillary fluid the matrix C is 
of the form 

/ 0* 0* \ 

C = -C* = -p e kk* 

\ Pe kk* J 

where p e is the equilibrium fluid density and is the zero- matrix 3x3. 

Here and later, for any vectors a, b we use the notation a*b for the scalar 
product (the line is multiplied by the column vector) and a b* for the tensor 
product (or a®b the column vector is multiplied by the line vector). Divergence 
of a linear transformation A is the covector div(A) such that, for any constant 
vector a, div(A) a = div(Aa). The identical transformation is denoted by /. 

In section 2 we present the multi-dimensional case in Eulerian coordinates for 
a particular form of the internal energy. In section 3 we consider in Lagrangian 
coordinates a one-dimensional case for the general form of internal energy. 



2 Governing equations in Eulerian coordinates 

The internal energy per unit volume of a capillary fluid is taken in the form 

i \ i \ c l w | 2 , ,\ 

e(p,»?,w) =e(p,fj)+-y- (4) 

where p is the fluid density, w = grad p (or w* = Vp), r] is the entropy per unit 
volume, c is the capillarity coefficient which is assumed to be a constant [Rocard, 
1952, Rowlinson and Widom, 1984]. The homogeneous energy e satisfies the 
Gibbs identity 

de = p dp + 9 dr) (5) 

where p = (e + P — Or))/p is the chemical potential, 6 is the temperature, 
de 

P = p— e is the thermodynamic pressure. By using Hamilton's principle, 

op 

governing equations of such a fluid were obtained by Casal (1972) (see also Casal 
and Gouin, 1985, Gavrilyuk and Shugrin, 1996). They are in the form 

p t + div j = 

* + div(5j)=o (6) 

j*+div ^+pJ + cww*j =0* 
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where j = pu, u is the velocity vector, index t is the partial derivative with 

Se ( ,. |w| 2 \ , S 

respect to time, p = p- e = P — c p divw H — , where — means the 

dp y 2 J dp 

variational derivative with respect to p. By using (4) and (5) we can obtain 
from (6) the energy conservation law 

( £+ 2^ + +diV (( £+ 2p +P ) u + c ( wdiv J "J div = 

(7) 

Since 

div |cww* - c ^pdiv w + I^j ^ 

= cw* div w + cw* (^—^j ~ cw * div w— cpV (div w) — cw* ^"^j—^ 
= — cpV (div w) 

the momentum equation reads 

j* t + div (0- +PI S j -cpV (div w) = 0* 

The gradient of the mass conservation law verifies another conservation law 

w* t + V (div j) = 
If we add an initial condition such that 

w*| t =o = Vp| t=0 

we can consider w as an independent variable. The fact that w* = Vp will be 
a consequence of the governing equations. Finally, we obtain equations (6) in 
the following equivalent non-divergence form 

p t + div j = 
Tit + div ^ j = 

or \ p , (8) 

j* t + div( — + PI\ -cpV(div w) =0* 
w* + V (div j) = 

The theory of capillary fluids is usually applied for van der Waals-like fluids. For 
such fluids the energy e (p, rj) is not convex for all values of p and ry. We suppose 
that we are in the vicinity of an equilibrium state (p e ,r] e ) where the energy 
function is locally convex. Let us introduce conjugate variables (q, 8, u, r) by 
the formula 

dE = d + |L + C ^L^J = ^ _ dp + dr] + u * d j+c w'dw (g) 

= qdp + 9drj + u*dj + r*dw 
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The Lagrange transformation of the total energy E is defined by 
IT = pq + r]6 + j*u + w*r - E = P + 

2c 

where the thermodynamic pressure P is considered as a function of q, 9 and u. 
Hence, in terms of the conjugate variables (q, 9, u, r) defined by equation (9), 
equations (8) can be rewritten in the following form 

dqj t \ dq J 

an\ / 9(nu) dndr\ () ( - w > 

du J t \ du dq dx J 

<m\ / o>(nu) 9nau\_ o 

dr J t \ dr dq dx J 

If the capillary coefficient c is zero, IT — P and we get the gas dynamics equation 
and the symmetric form of Godunov (1961). 

Multiplying equations (10) by q, 9, u and r, summing up all of them and 
using the identity 

rot (a x b) = [a, b] + a div b b div a 

where 

r ki da u dh 

[a ' b] = ^ b 

denotes the Poisson bracket, we get the conservation of the energy (7) in the 
form 

on e dn dn ^ an^ 

^ dq d9 du dr 

,. / <9 0lu) „d(Ilu) o>(nu) d(Uu) „ 
+ div q^- 1 + e^r 1 + u+— — ^ r-nu 
\ dq d0 du dr 

dU du^ dn dr\ _ 

dq dx dq dx J 

The system (10) admits constant solutions (p e , r/ e , u e , w e = 0). Since the gov- 
erning equations are invariant under Galilean transformation, we can assume 
that u e = 0. If we look for the solution of the linearized system proportional to 
g j(k x-At)^ we e q ua y on (3) 5 m w hich we have put 

*-■§-{{¥ 

dv \\dv 

C = -C* = 

and v* = (q, 9, u*, r*). Hence, eigenvalues A are real if A is positive definite. 
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3 One-dimensional barotropic case 

In mass Lagrangian coordinates (t, z) the governing equations are (see Gavrilyuk 
and Serre, 1995) 



v t -u z = 0, u t +p z = 0, p=- 



Se 


(de 






Sv 




dz 


\dv z )) 



, e = e{v,v z ) 



where v = — denotes the specific volume. This case is general: we do not suppose 
P 

a particular form (4) of the energy e. Consider an augmented system 



v t -u z =0 
w t -u zz = 
. de d ( de \ . 



dv dz \dw 
Let us define tt and the conjugate variables (a, r) as 

de de 
7r = — v + —— w = <jv + rw — e 
ov ow 

In terms of tt and (a, r, u) the system (11) reads 

(l)r"" =0 

1H- {(J- r z ) z = 



(11) 



In matrix form we get 
A 

where 



A = 



/ 


d 2 7T 


d 2 7T 




da 2 


dadr 




d 2 TT 


d 2 7T 




dadr 


dr 2 
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(13) 



5 



Equations (12) and (13) imply a dispersion relation of type (3), if we put B = 
kB\, C — k 2 C\. We note also that the system admits the energy conservation 
law 

where [r, u] = ur z — ru z . 

Remark. Analogous symmetric forms may be obtained for bubbly liquids, 
where the internal energy is a function not only of the density but also of the 
total derivative of the density with respect to time. 
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